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ABSTRACT 

The  Strategic  Defense  Initiative  has  generated  new 
interest  in  the  development  of  more  stable  space  structures. 
This  interest  has  increased  the  need  for  more  detailed 
knowledge  of  the  behavior  of  engineering  structures  under 
dynamic  loading.  Interests  lie  in  decreasing  the  amount  of 
vibration  by  both  passively  and  actively  damping  the 
structure.  M  means  exists  to  passively  damp  structures  by 
friction  damping  resulting  from  relative  slip  between  joint 
interfaces.  It  may  be  feasible  to  increase  the  damping  in  a 
structure  by  allowing  more  friction  damping  than  is  normal 
and  thereby  controlling  the  vibration  response. 

This  study  incorporates  friction  damping  in  a  one- 
dimensional  model.  Finite  element  techniques  are  used  to 
accomplish  the  numerical  analysis.  M  cl amped-cl amped  beam 
IS  used  as  the  physical  model.  The  mid-point  of  the  two 
element  beam  is  allowed  to  slip  in  rotation,  but  not  in 
translation.  because  the  one-dimensi onal  program  cannot 
handle  rotations  at  continous  nodes,  the  beam  is  modeled  by 
symmetry  about  the  joint  and  a  cantilever  beam  with  an 
applied  end  moment  is  studied. 

Results  for  the  response  of  a  beam  in  vibration  are 
presented  showing  displacement  of  the  joint,  relative 
rotation  at  the  joint,  and  relative  angular  velocity  at  the 
joint!  versus  time.  Various  clamping  pressures  and 


initiAi  loads  are  explored.  Diagrams  of  the  beam  shape  vs. 
time  shoM  the  shape  the  beam  takes  on  when  slip  occurs  at 
the  joint.  f-requency  calculations  show  that  the  period  of 
the  response  is  affected  by  clamping  pressure,  but  not  by 
the  initial  loading.  bnergy  loss  calculations  are  presented 


4 


for  various  clamping  pressures. 


CHAPTER  1 


INTRODUCTION 


1. 1  BACKORQUND 

There  has  been  an  .ncrease  in  interest  in  the 
stability  o-f  large  space  structures  due  to  the  Strategic 
Defense  Initiative  (SDI).  A  desire  exists  to  mount 
sensitive  equipment  (tracking  and  reflecting  devices) 
externally  on  these  structures.  The  rotation  of  these 
structures  could  cause  vibrations  to  occur  which,  if 
unchecked,  could  significantly  affect  the  accuracy  of  the 
mounted  equipment. 

Space  structures  operate  in  an  environment  offering  no 
aerod^naAiic  damping.  In  addition,  the  use  of  low  mass  and 
all  -welded  construction  methods,  lacking  sufficient  inherent 
damping,  decrease  the  ability  of  the  structures  to  reduce 
unwanted  vibrations  (19).  Thus,  expensive  and  complex 
damping  systems  _>  e  needed  to  achieve  the  desired  levels  of 
damping.  However,  friction  damping  is  a  method  available  to 
passively  damp  vibrations  in  structures,  thereby  reducing 
the  need  for  these  active  damping  systems. 

Friction  damping,  the  interfacial  slip  in  the  Joints 
of  a  structure,  is  the  major  contributor  to  the  inherent 
damping  of  a  structure.  Usually  over  90  percent  of  the 


damping  in  structures  takes  place  in  the  joints  (2).  Threi 


primary  reasons  are  given  for  not  relying  on  this  mechanism 
in  the  past.  First,  the  small  movements  between  the 
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surfaces  of  the  joint  can  cause  fretting  corrosion,  leading 
to  fatigue  and  possible  failure  of  the  joint  or  structure. 
Second,  by  decreasing  the  stiffness  in  the  joint  in  order  to 
permit  slip,  a  relative  loss  of  static  stiffness  of  the 
structure  is  realized.  Finally,  and  probably  the  most 
influential  reason,  friction  forces  in  the  joint  are  non¬ 
linear.  It  IS,  therefore,  difficult  to  determine  the  force 
transmitted  across  the  joint  and  to  predict  the  effect  on 
the  vibration  response  of  a  structure.  As  a  result  of  these 
disadvantages  of  friction  damping,  joints  are  normally 
clamped  tightly  to  prevent  interfacial  slip.  This  decreases 
the  damping  by  the  joints  to  a  minimum  (6). 

The  benefits  gained  by  friction  damping  in  joints  can 
no  longer  be  overlooked.  The  disadvantages  can  be  overcome 
by  improving  the  surface  finish  and  applying  joint  damping 
at  selected  joints,  as  opposed  to  all  the  joints.  Careful 
selection  of  joints  where  friction  damping  is  obtained  by 
slip  in  rotation,  but  not  in  translation,  need  sacrifice 
little  static  stiffness  (4).  Furthermore,  conservative 
results  may  be  obtained  by  linear  analysis  (5).  Therefore, 
it  may  be  feasible  to  increase  the  inherent  damping  in  a 
structure  by  allowing  more  friction  damping  than  is  pres¬ 
ently  allowed,  thereby  controlling  the  vibration  response. 
This  technique  can  help  to  avoid  complex,  active  damping 
systems.  Borne  of  the  advantages  would  be  lower  cost, 
achievement  of  high  damping,  the  ability  to  move  resonant 
frequencies,  and  the  ability  to  use  existing  joints  (5). 


V-  .• 


2 


1.2 


PURPOSE 


The  purpose  o-f  this  thesis  was  to  investigate 
Trictional  damping  e'f-fects  in  a  beam  with  a  rotational 
joint.  The  major  thrust  of  this  effort  was  to  take  an 
existing  finite  element  computer  program,  FEMID  (15),  and 
apply  friction  damping  to  change  the  vibration  response  of 
the  physical  model  and  thereby  show  the  capability  to 
control  the  vibration  response  and  resonant  frequency. 


1.3  General  ftpproach  and  Assumptions 

The  structure  that  presented  itself  most  readily  to 
serve  as  a  physical  model  for  this  study  was  a  two  element 
beam  clamped  at  both  ends  with  a  joint  located  at  the  middle 
node  of  the  beam.  A  sketch  of  this  set-up  is  shown  in  Fig. 
1.1a.  In  the  static  case,  when  the  beam  was  subjected  to  a 
vertical  load  at  the  midpoint,  the  beam  would  deflect  and  a 
moment  would  build  at  the  joint  due  to  a  frictional  moment 
opposing  slip  and  forcing  the  slope  to  remain  fixed  at  zero. 
When  the  moment  reached  the  value  of  the  friction  moment. 


any  further  increase  of  load  would  cause  the  joint  to  slip 
and  an  energy  loss  would  be  realized  due  to  the  relative 


rotation  o+  tne  two  interfaces.  H  new  displacement  and 
slope  would  b*?  obtained,  different  from  those  of  the  totally 
locked  joint  and  the  value  of  the  moment  at  the  joint  would 
equal  the  friction  moment  value.  t-igs.  1.1b  and  c,  show 
respectively  the  cases  where  the  moment  is  less  than  the 
friction  moment  and  when  the  beam  slips  because  additional 
load  IS  applied  once  the  friction  moment  is  reached.  In  the 
dynamic  case,  the  loss  of  energy  due  to  slip  would  be  seen 
as  a  decrease  in  the  amplitude  of  the  vibrating  beam  from 
cycle  to  cycle. 

A  literature  search  showed  that  interest  in  the  field 
of  friction  damping  in  joints  has  been  around  for  soma  time. 
However,  there  is  very  limited  vibration  data  available  in 
which  friction  forces  play  a  significant  role  (19).  The 
number  of  variables  associated  with  friction  damping  made  it 
difficult  to  predict,  so  it  was  avoided.  But  because  of 
bui ,  the  number  of  studies  undertaken  to  predict  the  effect 
of  friction  damping  in  joints  should  increase  significantly. 
Ihrough  the  judicious  use  of  assumptions,  the  scope  of  the 
problem  can  be  narrowed  so  that  a  starting  point  can  be 
obtained. 

Ihe  ability  to  predict  or  control  the  damping  produced 
by  partial  slip  between  surfaces  can  be  accomplished  in  the 
laboratory  on  only  very  simple  systems  ll4).  f-or  this 
reason,  macroslip  was  assumed  for  the  joiht  under  investi¬ 
gation.  In  other  words,  when  the  joint  was  slipping,  the 
total  contact  area  would  be  slipping.  This  leads  to  the 


figure  1.1a  TWO  ELEMENT  BEAM  MODEL 


FIGURE  1.1b  BEAM  WITH  M  <  Mqross 


FIGURE  1.1c  SLIPPED  BEAM 


next  assumption  regarding  the  clamping  pressure. 

The  clamping  pressure  was  applied  such  that  it  was 
evenly  distributed  across  the  contact  area.  I-f  a  bolt  is 
used  to  apply  the  clamping  pressure,  then  the  pressure  is  at 
a  maximum  near  the  bolt  location  and  decreases  as  you  go 
away  -from  it.  The  area  near  the  bolt  would  be  more  rigid 
and  would  resist  slip  more  than  the  outer  area.  The 
increased  pressure  near  the  bolt  would  mean  that  the 
friction  moment  would  be  greater  near  the  bolt.  In  order  to 
have  a  constant  friction  moment  across  the  contact  surface, 
which  the  macroslip  assumption  attempted,  a  uniform  clamping 
pressure  '*s  applied.  Another  asumption  was  required  to 
insure  a  constant  friction  moment. 

The  friction  force  applied  to  the  problem  is  propor — 
tional  to  the  clamping  pressure  through  the  coefficient  of 
friction,  fx.  Since  the  clamping  pressure  is  applied 
uniformly  over  the  contact  area,  the  coefficient  of  friction 
must  be  constant  to  provide  a  constant  friction  force  over 
the  entire  area.  This  assumption  of  constants  also  carries 
over  to  the  time-dependent  problem. 

Since  the  solution  of  a  vibrating  beam  is  a  time- 
dependent  problem,  the  clamping  pressure  and  coefficient  of 
friction  must  be  constant  in  time  or  the  friction  force  will 
vary.  The  first  step  is  to  find  a  solution  for  the  constant 
friction  force  or  too  many  variables  will  be  introduced  to 
solve  the  problem.  Hence,  the  friction  force  is  assumed  to 
be  constant  with  time.  This  leads  us  to  the  physical  beam 


model . 


M  simple  beam  model  that  is  available  is  the  buler 
beam.  I  he  theory  o-f-  the  Euler  beam  states  that  a  plane 
cross-sectioh  drawn  perpendicular  to  the  cehterline  o+  ah 
uhbeht  beam  will  remaih  perpendicular  to  the  cehterline 
after  the  beam  behds  (V)  .  Ihe  limoshehko  beam  does  not 
require  the  cross-secti oh  to  remaih  perpendicular  to  the 
beam  cehterline  and  so  more  closely  approximates  the  exact 
beam  <V) .  The  dimensions  of  the  clamped-clamped  beam  were 
chosen  such  that  the  results  obtaihed  for  the  first  mode 
would  closely  approximate  that  of  the  limoshehko  beam. 
Lohsequently ,  the  simple  thin-beam  Euler  model  will  give 
adequate  results  for  the  first  mods. 

Einally,  a  modification  had  to  be  mads  to  the  clamped- 
clamped  beam  model.  A  one-dimensional  finite  element  pro¬ 
gram  was  used  to  model  the  physical  beam  and  provide  the 
numerical  analysis  for  the  dynamics  of  the  beam.  The  one- 
dimensional  program  did  not  allow  for  independent  rotations 
of  the  elements  at  continuous  nodes|  it  forced  the  elope  at 
the  middle  node  to  remain  zero.  but  for  damping  to  take 
place,  relative  rotations  had  to  occur  at  the  node  to  dissi¬ 
pate  energy.  bo  the  beam  was  changed  by  symmetry  about  the 
joint  to  a  cantilever  beam  with  an  applied  end  moment.  The 
slope  would  be  fixed  whenever  the  moment  was  lass  than  the 
frictional  moment,  and  would  be  allowed  to  change  when  the 
value  of  the  moment  reached  the  friction  moment  value.  Ihis 
allowed  the  rotations  required  to  calculate  energy  loss. 


CHAHIth  2 


IHEOHY 

2. 1  Loul omb  Uampi nq  T heory 

CJoiilomb,  or  dry  -friction  damping  arises  Mhenever  two 
bodies  are  allowed  to  slide  or  rub  against  one  another.  For 
any  sliding  to  take  place  there  must  be  a  -force  acting  which 
overcomes  the  resistance  caused  by  friction.  The  friction 
force  IS  parallel  to  the  surface  and  proportional  to  the 
force  normal  to  the  surface.  The  classical  law  of  sliding 
friction  states  that  this  frictional  force  is  independent  of 
the  contact  area  and  the  magnitude  of  the  velocity,  as  long 
as  sliding  exists  (9).  The  friction  force  opposes  the 
relative  motion,  and  thereby  continuously  absorbs  energy  so 
long  as  the  relative  velocity  exists.  The  force  of  friction 
is  of  constant  magnitude  and  as  long  as  the  forces  acting 
(namely,  the  inertia  force  and  the  restoring  force)  are 
sufficient  to  overcome  the  friction  force,  damping  will 
continue.  When  the  forces  become  too  small,  the  damping 
stops  as  sliding  ceases. 

This  can  be  illustrated  by  the  motion  of  the  simple 
spring-mass  system  of  Fig.  2.1. 

The  mass  is  set  into  motion  by  pulling  it  in  the 
positive  direction  and  releasing  it.  The  friction  force 
acts  in  the  opposite  direction  of  the  velocity  of  the  mass. 
It  is  proportional  to  the  force  acting  normal  to  the  contact 


F=-mW  =  Mhcy 

FIGURE  2.1  SIMPLE  SPRING-MASS  SYSTEM 

sur+ace,  in  this  case  the  ^aight  o-f  the  mass.  The  constant 
of  proportional  ity  is  the  coefficient  of  friction,  ft.  This 
constant  depends  only  on  the  roughness  of  the  sliding 
surfaces.  The  force  of  friction  is  therefore  represented  by 

Fd  «  *  pmg  (2-1) 

The  equation  of  motion  for  this  system  can  be  written 
as 

mx  +  Fcisgn(x)  +  kx  =  O  (2-2) 

where  the  symbol  "sgn"  represents  a  function  having  a  value 
of  +1  if  X  IS  positive  and  -1  if  x  is  negative.  The 
equation  of  motion  for  this  system  is  non-linear,  but  it  can 
be  separated  into  two  linear  equations,  one  for  the  positive 
velocity  and  one  for  the  negative  velocity  (12). 


•• 

mx 

+  kx  =  -Fc 

* 

X 

> 

0 

(2-3a) 

•• 

mx 

+  kx  »  Ft, 

• 

X 

< 

0 

(2-3b) 

The  non-linearity  consists  of  switching  between  the  two 
linear  equations.  The  switching  does  not  occur  as  an 


explicit  function  of  time,  but  is  determined  by  the 


response.  Ihere+ore,  it  could  occur  at  various  times, 
depending  on  the  torcing  'function  and  the  initial 
conditions.  This  means  the  equation  of  motion  must  be 
solved  for  one  time  interval  at  a  time,  depending  on  the 
sign  of  X. 

If  the  mass  is  displaced  initially  in  the  positive 
directir'.,  the  velocit/  Mill  be  negative  at  first,  so  the 
first  equation  is  considered.  It  can  be  rewritten  in  the 
form 

•*  2  2 

X  +  Wr,x  «  w,,fd  (2-4) 

k 

where  w.,  *  -  (Natural  Frequency  of  Vibration) 

m 

F<. 

fd  *  - 
k 

The  initial  conditions  are  x(0)">Xc  and  k(0)>bO  bo  the 
solution  is  (12) 

x(t)  ■  (Xo  -  fd)co5(Wr.t)  +  fd  (2-5) 

This  represents  harmonic  oscillation  and  is  valid  for 
O  <  t  £  ti,  where  ti  is  the  time  at  which  the  velocity 
reduces  to  zero  and  the  motion  is  about  to  reverse.  To  find 
the  velocity  the  equation  is  differentiated  with  respect  to 
time  to  obtain 

x(t)  ■  -Wr,(Xo  -  f <j) sin (Mr.t )  (2-6) 

BO  that  the  lowest  non-trivial  solution  satisfying  the 
initial  condition  is  ti*  s/Wr..  If  x(ti)  is  large  enough  to 
overcome  the  static  friction,  the  mass  will  have  a  positive 
velocity  and  then  must  satisfy  the  other  linear  equation  of 


i? 


motion.  Thus  the  following  equation  must  be  solved 

X  +  Wr,X  ■  -  W.^-fd  (2-/) 

with  the  initial  conditions  of  x(ti)*  -(x^  -  2fd)  and 
x(tilsO.  fhe  solution  is  (12) 

x(t)  ■  (Xd  -SfdlcoBlWdt)  -  +ci  (2-a) 

This  solution  IS  valid  for  the  values  of  t»  i  t  s  ts  where 
tz  is  the  value  at  which  the  velocity  again  goes  to  zero. 

The  value  of  ts  is  found  to  be  2«/Wd.  The  procedure  is  then 
repeated  for  t  >  tz,  until  the  motion  stops.  However,  a 
pattern  is  seen  to  emerge. 

Over  each  half-cycle  the  motion  consists  of  a  constant 
component  and  a  harmonic  component  with  the  frequency  equal 
to  the  natural  frequency  of  the  spring-mass  system,  where 
the  duration  of  every  half-cycle  is  equal  «/Wr,.  The  average 
value  of  the  solutions  switches  between  fd  and  -fd  end  at 
the  end  of  each  half -cycle  the  displacement  is  reduced  by 
2fd  >■  2f<mg/k.  This  leads  to  the  conclusion  that  in  coulomb 
damping  the  amplitude  decays  linearly  with  time,  not 
exponentially  as  in  viscous  damping  (12).  A  plot  showing 
the  decay  as  a  function  of  time  for  this  simple  system  is 
seen  in  Fig.  2.2.  Fig.  2.2  also  shows  that  the  curve 
oscillates  about  the  values  of  fd  and  -fd  and  if  the 
restoring  force  falls  below  this  value,  the  block  will  stop 
sliding  at  some  position  on  either  side  of  zero,  depending 
on  the  initial  displacement. 

One  may  associate  coulomb  friction  with  rotational 
slip  by  looking  at  the  force  of  sliding  friction  coming  from 
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FIGURE  2.2  LINEAR  FRICTION  DAMPING  DECAY  (12) 


a  -frictional  shear  stress  at  the  contact  area  that  is  equal 
to  u  times  the  normal  pressure.  This  gives  an  equivalent 
result  to  the  simpler  cases,  such  as  the  sliding  block 
previously  mentioned.  This  case  aids  in  the  solution  o-f 
more  complicated  systems,  such  as  curved  cohtact  sur-faces  or 
non-uni-form  pressure  or  velocity  di str i buti ohs  (10). 

To  show  this  case,  assume  a  frictional  moment  arises 
from  the  flat  end  of  a  circular  rotating  shaft  of  radius  R 
being  pressed  against  a  plane  surface  with  a  total  force  N 
(lU).  Also  assuming  a  uniform  normal  pressure  or  compres¬ 
sive  stress  on  the  contact  area  A,  we  obtain  a  uniform 
frictional  stress  of  magnitude  (lo) 

N  ^N 

7^,  =  ^  « -  (2-9) 

A  <«R- 

which  IS  everywhere  in  a  direction  normal  to  a  radial  line 
drawn  from  the  center  of  the  circular  contact  area.  From 
Fig.  2.3  it  IS  seen  that  the  moment  due  to  an  annular 


-  w  y  •  V  -^5  ■>  ->  "k  •  U  *1  •  '. 


element  o-f  width  dr  and  area  2'Krdr  is 


»  i 


pN 

dM  =  7^2'Kr^dr  =  —  2'Kr==dr 

ifR^ 


(2-10) 
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If 

p  =  0,  the 

contact 

surface 

IS 

said  to 

be  perfectly 

m 

smooth.  If,  on  the  other  hand,  »<  =  ••»  a  perfectly  rough 
surface  exists.  The  first  case  corresponds  to  the  no  clamp 
case  in  which  the  friction  force  is  zero.  The  second  case 
corresponds  to  the  relative  velocity  at  the  point  of  contact 
being  zero  and  is  referred  to  as  the  no  slip  case  (10). 

As  an  extension  to  this  case,  Richardson  and  Nolle 
assumed  that  an  external  moment  applied  to  a  rotary  joint  is 
transferred  between  the  two  surfaces  by  shear  due  to 
friction  (16).  Fig  2.4  depicts  the  transfer  at  the  friction 
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The  moment  is  applied  -from  zero  to  a  value  o-f  MriAx . 

The  inner  bounding  radius  can  be  calculated  from  equilibrium 
of  moments. 

M  ■  \  2‘Kr^pPdr  yields  a^  “  R-^  -  3M/2n#4p  (2-12) 


where  P  “  Clamping  Pressure 

For  gross  eliding,  the  slipped  region  covers  the 
entire  contact  area  (a  “  0) ,  so  the  moment  for  gross  slip  to 
occur  IS 

MoFtosB  2s^PR^/3  (2—13) 

For  any  applied  moment  below  McaFtoss  the  relative  slip  is  in 
the  microslip  stage.  At  lioFtass,  macroslip  would  be  said  to 
to  occur.  The  value  of  Mcsbosb  can  be  varied  as  a  function 
of  clamping  pressure.  The  low  end  of  the  macroslip  region 
occurs  at  low  clamping  forces  and  is  known  as  the  no-clamp 
case.  As  the  moment  increases  to  Moross.  the  amount  of 
damping  increases  until  Moross  is  reached.  Above  Morobs, 
damping  decreases  in  the  microslip  region  until  the  clamping 
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FIGURE  2.6  DAMPING  VS  CLAMPING  PRESSURE  <11 ) 


Timoshenko  beam  does  not  require  this  to  happen  and 
therefore  models  the  real  beam  better.  In  an  effort  to 
estimate  the  amount  of  error  thin  beam  theory  would  enter 
into  the  problem,  it  was  compared  to  Timoshenko  theory. 

Fig.  2.7  compares  a  non-dimensional  frequency  for  the 
Timoshenko  beam  vs  the  classical  beam  (thin  beam).  The 
length  to  height  of  the  physical  beam  model  gives  an  L/r 
value  of  68.  This  compares  quite  well  with  the  Timoshenko 
beam  for  the  first  mode,  as  shown  by  the  *  in  Fig.  2.7. 
Errors  may  show-up  for  higher  vibration  modes,  but  the  thin- 
beam  theory  will  be  adequate  for  predicting  the  first 
vibration  mode. 

2.3  Fini te  Element  Theory  for  T i me-Dependent  Problems 

In  time-dependent  problems,  the  undetermined 
parameters  are  assumed  to  be  functions  of  time  and  the 
approximation  functions  are  assumed  to  be  functions  of 
spatial  coordinates.  Two  stages  of  solution  must  be 
undertaken,  both  involving  approximate  methods.  First,  the 
spatial  approximation  must  be  considered  and  then  the  time 
approximation.  This  procedure  is  known  as  semidiscrete 
approximation  (in  space)  (IS).  The  spatial  approximation 
leads  to  a  set  of  ordinary  differential  equations  in  time. 

In  structural  dynamics,  these  equations  quite  often  contain 
second-order  time  derivatives  of  the  dependent  variables. 

The  spatial  approximation  forms  the  equations  of  motion  in 


matrix  form 


CM3fVj  +  CK3tx>  «  {F>  0  <  t  <  to  (2-14) 

where  LM3  “  mass  matrix 

CK]  B  stiffness  martix 
<x>  =  displacement  vector 
{F>  s:  applied  loads  vector 

This  equation  must  be  further  approximated  to  get  a  set  of 
algebraic  equations  which  can  then  be  solved.  There  are 
several  schemes  available  for  solving  the  above  equation. 

The  most  common  one  is  the  Newmark-Beta  direct  integration 
method  (15).  The  Newmark-Beta  method  approximates  the  first 
time  derivative  {x>  and  the  function  {x>,  which  is  time 
dependent,  at  the  (n  -*■  Dst  time  step  (at=  constant)  as (15) 

<X>r.^l  ■  <X>r.  +  C  (  l-T)  {x'lr,  +  TtVloH-llAt  (2-15a) 

tx>n+l  -  tx>o  +  <x>oAt  +  C(‘4  -  B)<x>o  +  BCx>r,^i3  (At)*(2-15b) 
wh  e  7  and  B  control  the  accuracy  and  stability  of  the 
scheme  and  subscript  n  shows  that  the  solution  is  computed 
at  time  t  ■  tr>.  The  values  of  t  *  4  and  B  ■  >4  are  normally 
chosen  because  they  provide  for  an  unconditionally  stable 
system.  This  means  that  the  size  of  the  time  step  is 
goverened  by  the  need  for  accuracy  rather  than  stability. 
This  scheme,  known  as  the  constant  average  acceleration 
method  (15),  is  very  valuable  in  computing  the  response  of 
systems  in  which  damping  is  introduced  externally  because 
the  method  itself  does  not  introduce  any  damping.  Equations 

(2-14)  and  (2-15)  can  be  combined  and  rearranged  to  get  (15) 

A  A 


(2-16) 


A 

where  tK3  ■  CKD  +  «oCM3 

CF>r.-n  =  CF>ri-.-i  +  LM3(ao'Cx>r.  +  ai(x>n  +a2Cx>r.) 

Given  initial  conditions  ai  <x}o,  <><>0,  and  <x*>o,  equation 
(2.16)  can  be  solved  repeatedly  for  fx>.  1-f  the  initial 
value  of  the  acceleration  is  not  given,  which  is  usually  the 
case,  equation  (2-14)  can  be  used  to  compute  <x*>  at  t  *  0. 
After  calculating  <xJ  at  tr.-^»  »  (n+Dat,  the  first  and 
second  derivatives  (velocity  and  acceleration)  of  fx>  at 
tr^^i  can  be  found  by  rearranging  equation  (2-15)  (IS) 


fX*>r,-.-l 

«  •» 

<x>r.  -  a^fxlr. 

(2-17a) 

• 

C  X  ^  ri  1 

*  <x>„ 

•• 

+  asfxlr^  +  a^fxlr. 

1 

(2- 17b) 

Where 

ao  = 

l/C6(At)=l  a. 

(1/26)  -  1 

as  = 

(1  -  7) At  a« 

»  'TAt 

Given  initial  conditions  of  Cx>o,  <x>o,  and  fx>o,  equation 
(2.16)  can  be  solved  repeatedly  for  fx>  and  its  time 
derivatives  at  any  time  t  >  O. 

Newmark-Geta  is  unconditionally  stable,  as  previously 
mentioned,  which  means  the  solution  does  not  grow  without 
limit.  However,  if  the  right  time  step  is  not  chosen  the 
results  may  be  inaccurate.  An  estimation  for  the  time  step 
is  given  by  (15) 

At  ■  (2-18) 

where  is  the  smallest  period  of  vibration  given  by  the 

problem.  This  assures  not  only  a  stable  solution  for  any  at 
but  also  accuracy  of  results.  The  smallest  period  in  the 
cl amped-cl amped  beam  occurs  when  the  joint  is  fully  clamped. 
The  period  was  obtained  from  the  unmodified  FEMID  program. 


The  value  was  0.00/4074  sec.  Substituting  this  number  into 
Eqn.  <2-l8> ,  a  at  oT  0.0023579  was  computed.  Therefore,  any 
value  of  at  below  this  number  would  guarantee  accuracy  of 
results.  The  value  of  at  =  0.0001  was  chosen  for 
calculation  in  the  Newmark  Beta  method  in  an  effort  to 
obtain  increased  accuracy. 

Fig.  2.8  highlights  the  main  steps  of  the  Newmark  Beta 
time  integration  method.  Fig.  2.9  compares  the  accuracy  of 
the  Newmark  Beta  method  with  other  numerical  analysis 
methods  and  the  exact  solution.  The  problem  being  solved  in 
Fig.  2.9  is  a  two  element  axial  rod  clamped  at  one  end.  The 
top  plot  shows  the  axial  displacement  of  the  middle  node  and 
the  bottom  plot  shows  the  axial  displacement  of  the  node  at 
the  free  end.  Notice  that  while  the  other  numerical  methods 
display  some  form  of  amplitude  decay,  the  Newmark  Beta 
method  does  not.  However,  Newmark  Beta  does  exhibit  some 
period  elongation,  and  this  should  be  considered  if  accurate 
frequency  calculations  are  required. 
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FIGURE  2.B  NEWMARK  BETA  INTEGRATION  METHOD 


A.  Initial  Lalculationa 

1.  Form  K,  and  M. 

2.  Initialize  and 

3.  Select  time  step  at,  i  and  B  (t  >  *4,  fit  >4(0.5  +  t)^), 
and  calculate  : 

ao  =  l/(BAt^)  ;  ai  «s  l/ifiAt)  ;  a^  *  (1/2B)-1 

as  =  At  ( 1 — i)  ;  a.*  “  TAt 

A 

4.  Form  the  matrix  K  “  K  +  aoM. 

B.  For  each  time  step 

1.  Calculate  the  cf-fective  loads  at  time  t  At  t 

F«;...a*:  “  F«i-»-At  +  MtaoXt  +  azXfc  +  asX^J 

2.  Solve  for  the  di spl acements  at  t  At  t 

A  A 

KX-t-.-At:  “  F-tL-.-Ati 

3.  Calculate  the  accelerations  and  velocities  at  t  At  : 

•  «  •  •• 

-  Xt>  -  *iKt:  -  azXii 

n  • 
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CHHHIbH  3 


THE  DttnUNb  MUUhL 

3. 1  FHb  BbHM  IN  blHI lU  BbNUlNb 

Ihe  physical  dimensions  oT  the  two  element  beam  arrangement 
can  be  seen  in  Fig.  3.1.  However,  because  oT  the 
limitations  of  the  FEMIU  computer  program,  the  actual  beam 
model  that  was  used  in  the  computer  analysis  was  a 
cantilever  beam  with  an  applied  end  moment.  The  beam  was 
made  of  steel;  the  material  properties  and  related 
calculations  to  be  used  in  the  computer  analysis  can  be 
found  in  Appendi>i  1. 


FIb'URE  3.1  PHYSICAL  BEAM  DIMENSIONS 
The  computer  program  allowed  certain  boundary 
conditions  to  be  fixed,  so  initially  it  was  assumed  that  no 
slip  would  take  place  and  the  slope  at  the  tip  was  fixed  at 


zero.  As  a  vertical  load  Mas  applied,  a  moment  was 
generated  to  keep  the  slope  fixed  at  zero.  It  the  value  o+ 
Mqrobs  was  not  equaled,  then  no  slip  took  place  and  the  beam 
would  behave  as  a  normal  cl amped-cl amped  beam  in  bending. 
When  the  value  ot  the  moment  reached  the  value  ot  MGRQSS, 
the  next  increment  ot  load  would  cause  the  joint  to  slip 
because,  due  to  friction,  the  joint  cannot  hold  any  more 
moment  than  the  value  o+  MtaRoss.  At  this  point,  the 
boundary  conditions  were  changed  and  the  restriction  that 
the  slope  remain  zero  was  removed.  A  moment  equal  to  the 
value  ot  McsBasa  was  applied  opposite  to  the  direction  ot 
impending  motion,  and  slip  occured  due  to  the  change  in 
slope.  A  new  value  ot  displacement  and  slope  were  found, 
different  than  the  value  tor  the  cl amped-cl amped  beam.  Fig. 
3.2  shows  the  shape  ot  the  beam  and  equation  that  is  solved 
tor  the  static  beam  where  the  moment  does  not  equal  Mcsbogs. 
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FIGURE  3.2  UNSLIKPEI)  STATIC  BEAM 


Unce  slip  occurred,  a  new  equation  was  solved  because  ot  the 
boundary  condition  change.  This  is  seen  in  Fig.  3.3. 
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FIGURE  3.3  SLIPPED  STATIC  BEAM 
The  solution  o-f  the  beam  in  either  case  can  be  verified  by 
superpositon  once  the  value  of  the  moment  is  known.  The 
equations  of  superposition  for  beam  bending  by  an  applied 
load  and  moment  at  the  tip  of  a  cantilever  beam  are  (7) 

PL-^  ML=  PL^  ML 

„  =  _  -  +  -  ANO  e  =  -  _  —  (3-1) 


PL-i 

ML  = 

PL^ 

ML 

-  +  - 

AND 

e  =  - 

_  — 

3EI 

2EI 

2E1 

El 

where  w  =  displacement 
0  =  slope 
P  =  applied  load 
M  =  applied  moment 
L  =  length 

E  =  modulus  of  elasticity 
I  =  moment  of  inertia 


-  -.  --J 
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3.2  THE  BEAM  IN  DYNAMIC  BENDING 

Once  the  static  solution  was  found,  the  displacement 
vector  was  used  as  initial  conditions  for  a  dynamic  problem. 
The  initial  boundary  conditions  are  such  that  the  slope  is 
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-fixed  at  the  value  calculated  -from  the  static  problem.  The 
beam  is  then  set  into  motion  by  removing  the  initial 
vertical  load.  The  beam  begins  to  move  with  some  velocity, 
W2  (see  Fig.  3.4).  Initially,  the  moment  that  was  held  at 
the  joint,  the  value  of  Mofioss,  is  relaxed  because  of  load 
removal.  However  because  the  slope  is  not  allowed  to  change 
the  moment  begins  to  build  again.  The  following 
differential  equation  now  applies 

Mx*  +  Kx  =  F  (3-2) 

where  H  «  global  mass  matrix 

K  «  global  stiffness  matrix 
F  s  applied  loads  (F  «  <  0,  0  >^> 

Fcemember  that  the  static  equation  was  just 

Kx  -  F  (3-3) 

The  motion  of  the  beam  requires  the  inclusion  of  the 
acceleration  terms  in  the  problem.  This  presents  two  new 
terms  in  the  calculation  of  the  moment  to  compare  against 
McaFtcise.  The  moment  calculated  was  therefore 

M  ■  H^sWz  ^  M*.*02  +  K.43W2  +  K««S2  (3“4) 

where  the  displacement  and  acceleration  vectors  are  obtained 
by  solving  Eqn. (3-2)  by  the  Newmark  Beta  method.  Appendix  2 
shows  how  the  boundary  conditions  are  applied  to  determine 
which  terms  will  appear  in  the  dynamic  moment  calculation. 
Now  when  this  moment  equals  the  value  of  MoFtoea,  slip  will 
occur  with  the  next  incremental  increase  in  the  moment.  It 
can  be  seen  that  initially  62  is  zero,  but  after  slip  occurs 
an  angular  acceleration  term  will  appear  and  figure  into  the 


moment  equation.  Ihe  n:?w  equation  that  must  be  solved  once 
slip  occurs  IS 

mV  +  Kx  =  K  (3-5) 

where  F  =  applied  loads  (F  =  1  O,  Mcsrobs  >'). 

Fig.  3.4  shows  the  conditions  that  exist  once  slip  occurs. 


h 
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FIl3UF<t  3.4  DYNAMIC  BEAM  WITH  8=  <  0 
□nee  the  slip  mechanism  is  turned  on,  some  method  must 
be  available  to  turn  it  oft.  In  friction  damping  the 
friction  force  acts  opposite  the  direction  of  motion  of  the 
system.  In  this  case  the  motion  is  rotational,  so  the 
friction  moment  opposes  the  angular  velocity  of  the  joint. 
Iherefore,  the  angular  velocity  should  play  some  part  in 
monitoring  when  slip  at  the  joint  had  stopped.  When  slip 
ceases,  the  angular  velocity  should  be  equal  to  zero.  At 
this  point,  the  slope  was  fixed  and  the  friction  moment 
removed.  The  beam  would  then  begin  to  move  in  the  opposite 
direction  and  the  process  would  begin  again.  Fig.  3.5  shows 
this  case.  The  moment  at  this  point  is  positive,  but  it 
decreases  with  time  until  it  becomes  negative.  When  the 
absolute  value  of  the  moment  on  this  stroke  equaled  the 
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VAIUS  O^f  MtTiRcjoo  5  the  slope  Mas  aiioMed  to  change  and  the 
friction  moment  Mas  applied  in  the  opposite  direction  of  the 
motion.  This  is  illustrated  by  Fig.  3.6.  Th«  procesa  WBuia 
continue  until  the  moment  no  longer  equalled  the  value  of 
Mcsrosb  and  the  beam  would  oscillate  from  then  on  with  some 
permanent  slope.  Fig.  3.7  shows  the  logic  implemented  to 
solve  for  the  damping  at  the  joint  of  a  vibrating  beam.  The 
major  modifications  were  made  in  the  STRESS  subroutine  of 
the  FEMID  program  (15).  The  actual  Fortran  code  for  this 
subroutine  can  be  seen  in  Appendix  3. 


FIGURE  3.5  DYNAMIC  BEAM  WITH  02  =  0 


FIGURE  3.6  DYNAMIC  BEAM  WITH  02  >  0 


f-ILiUKt  ^./  LUUIL  Uh  I  Ht  t-KiCliUN  DAMPING  MhUHANISM 


hor  the  initiei  loading  o+  f-ig.  1.1c,  the  flow  of  the 
damping  mechaniEm  would  be: 

A.  hor  each  time  step  when  not  slipping 

1.  UalCUlate  Mean  OSB 

McaRoss  =  / S  “  Constant 

u:.  calculate  n 

•  •  ** 

M  *  k^3W3  + 

si.  Compare  M  vs  Moross 

a.  If  |n|  ^  Mqrosb  (so  on  to  next  time  step 

b .  If  1^1  ^  Mqrobs 

1.  Kemove  fixed  slope  condition 

2.  Apply  tip  moment  equal  to: 

a.  MtaROSES  If  n  ^  U 

b .  “McaROSB  1  f  M  <■  U 

3.  (so  on  to  next  time  step  in  part  B. 

b.  Por  each  time  step  that  slip  is  occurring 

1  •  I  M  j  ~  Moross 

2.  Check  angular  velocity,  Bz 
a.  If  n  >  0 

1.  If  8z  <  <J  (io  on  to  next  time  step 

2.  If  02  >  O 

a.  Fix  slope  at  present  value  (62  >«  Constant) 

b.  Remove  tip  moment 

c.  bo  on  to  next  time  step  in  part  a. 
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CHAPTER  4 


RESutjrs 

4 .  1  TH[E  STATIC  BEAM 

The  response  of  the  static  beam  in  bending  was  readily 
calculated  once  the  friction  moment  Mas  obtained.  Ihe 
method  consisted  of  applying  superposition  when  the  value  of 
Mas  reached.  The  applied  load  and  moment  Mould  then 
give  a  displacement  and  slope  at  the  end  node.  The  values 
for  the  the  static  problem  Mere  used  as  the  initial  dis¬ 
placement  vector  for  the  dynamic  problem.  Since  this  is 
Mhere  the  interest  lies,  emphasis  is  placed  on  the  results 
of  the  vibrating  beam. 

4.2  THE  DYNAMIC  BEAM 

The  beam  Mas  set  into  motion  Mith  the  slope  at  the 
joint  fixed  by  removing  the  initial  load.  When  the  dynamic 
moment  reached  the  value  of  Mcbriobb,  the  boundary  condition 
that  kept  the  slope  fixed  Mas  removed  and  a  friction  moment 
equal  to  the  value  of  Morobb  Mas  applied  in  opposition  to 
the  value  of  the  angular  velocity.  The  joint  Mas  alloMed  to 
slip  until  the  angular  velocity  reached  zero.  Then  the 
slope  Mas  fixed  at  the  value  it  had  reached  at  that  time  and 
the  moment  Mas  removed.  This  process  continued  until  the 
slope  remained  fixed  continuously  and  no  more  slip  took 
place.  The  values  that  Mere  deemed  important  Mere 
displacement,  slope,  and  angular  velocity.  All  three  of 
these  quantities  Mere  plotted  versus  time  to  determine  hOM 
the  beam  modal  Mas  responding  as  a  function  of  time.  Figs. 


4.1  through  4.4  •re  plots  o+  the  translational  displacements 
versus  time  tor  various  clamping  pressures.  Note  how  the 
maximum  displacement  decreases  as  the  clamping  pressure 
increases  tor  an  initial  applied  load  ot  1U,UUU  N.  The 
maximum  displacement  occurs  at  zero  clamping  pressure  and 
the  minimum  at  the  tully  clamped  case.  At  low  clamping 
pressures,  damping  takes  place  over  a  longer  time  period 
because  ot  the  small  value  ot  Morobs,  while  at  the  higher 
clamping  pressures,  damping  takes  place  tor  only  one  or  two 
halt-cycles. 

Figs.  4.5  and  4.6  show  plots  ot  slope  at  the  joint 
versus  time  tor  two  clamping  pressures.  At  the  lower 
clamping  pressure,  the  slope  was  seen  to  level  ott 
periodically  and  then  continue  after  it  had  been  fixed  for  a 
few  time  steps.  This  was  confusing  until  the  angular 
velocity  versus  time  plots  were  examined.  These  plots  are 
shown  in  Figs.  4.7  to  4.9. 

The  angular  velocity  versus  time  plot  for  zero 
clamping  pressure.  Fig  4.7,  exhibited  an  unexpected 
response.  Note  that  a  high  frequency  response  is  superim¬ 
posed  on  a  lower  frequency.  It  was  determined  that  two 
modes  ot  vibration  were  showing  on  this  plot.  The  lower 
frequency  trace  was  the  first  mode  of  the  beam  and  the 
higher  frequency  in  between  was  the  second  mode  due  to  the 
rotation  of  the  end  of  the  beam.  Whan  the  clamping  presura 
was  increased  to  b.0k06  N/m^,  both  modes  ware  seen  to  damp 
out,  as  would  be  expected.  The  higher  frequency  second  mode 


IRE  4.6 


FIGURE  4.7  Angular  Velocity 


FIGURE  4.8  Angular 


so.  O't 


FIGURE  4.9  Afigular  Velocity  vs.  Time 


damped  out  fir«t,  and  the  tiret  mode  then  continued  on 
because  there  was  no  structural  or  viscous  damping  modelled. 

Frequency  was  plotted  as  a  function  of  clamping 
pressure  and  the  lower  clamping  pressures  gave  lower 
frequencies.  As  the  clamping  pressure  was  increased,  the 
beam  became  more  stiff,  and  so  the  frequency  increased. 

This  IS  seen  in  Fig.  4.10.  Note  that  the  range  of 
frequencies  is  bracketted  by  the  values  obtained  for  a 
cantilevered  beam  at  the  low  end  and  a  fully  clamped  beam  at 
the  high  end. 

The  slope  at  which  the  joint  finally  locks  up,  the 
permanent  set,  is  a  function  of  the  initial  load  and  the 
clamping  pressure.  It  can  lock  up  at  a  positive  or  negative 
slope.  Figs.  4.11  through  4.13  depict  Wz,  8s,  and  63  versus 
time  respectively  for  an  initial  load  of  10,000  N  and  a 
clamping  pressure  of  V.SEOa  N/m-^.  The  permanent  set,  which 
was  a  negative  slope  for  most  of  the  previous  runs,  is  now 
set  at  a  positve  slope.  The  same  effect  could  have  been 
achieved  for  a  constant  clamping  pressure  and  varying 
initial  load.  The  beam  then  oscillates  about  some 
equilibrium  point  other  than  zero.  The  equilibrium  point 
could  also  be  either  positive  or  negative,  again  based  on 
the  initial  conditions. 

The  shape  of  the  beam  was  plotted  as  a  function  of 
time.  Figs.  4.14  and  4.15  show  the  shape  of  a  fully  clamped 
joint,  BO  that  the  slope  at  the  joint  remains  zero.  Figs. 


4.16  thru  4.18  are  for  a  clamping  pressure  of  B.0E06  N/m^. 
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FIGURE  4.13  Angular  Velocity 


Fully  Clammed  Beam 
Initial  Load  is  Id, 000  IJ 
Timestep  is  0.0001  sec 


Clar.^iiic  .-rfe^Eure  ic.  8.U£,Jt  n/. 
Initial  Load  is  10,000  N 
Timestei  is  0.0001  sec 


FIGURE  4.18  Beam  Shape  vs.  Time 
Third  Half-Cycle 


The  slope  and  displacement  are  seen  to  vary  si gni -f i cantl y 
because  the  small  value  o-f  Moross  is  reached  almost 
immediately.  Figs.  4.1V  and  4.2u  are  for  a  clamping 
pressure  of  B.OfciOU  N/m^.  The  displacement  is  damped  a  large 
amount  in  the  first  half-cycle,  and  the  remainder  is  damped 
in  the  next  half-cycle.  The  beam  then  vibrates  with  a  fixed 
slope  of  -u.UlOlti  rads  at  the  joint. 

Fig.  4.21  is  a  plot  of  moment  versus  slope  for  a 
moderately  large  clamping  pressure  of  B.OEiOB  N/m^  and  shows 
a  hysteresis  loop  for  the  damping  taking  place  at  the  joint. 
If  the  joint  were  excited  by  a  cyclic  load,  the  loop  would 
be  closed.  Damping  takes  place  on  the  horizontal  lines  of 
the  loop,  where  the  moment  is  constant  and  the  slope  is 
changing.  This  allows  the  energy  loss  to  be  calculated  by 
multiplying  the  moment  by  the  change  in  slope.  Remember 
that  before  damping  begins,  the  dynamic  moment  must  equal 
the  gross  moment.  The  program  was  unable  to  reach  the  gross 
moment  exactly,  so  some  overshoot  of  the  gross  moment  is 
seen.  These  are  the  spikes  just  before  the  line  turns 
horizontal.  One  other  spike  is  seen  on  the  upper  horizontal 
line.  This  is  a  point  where  the  angular  velocity  changed 
sign  and  so  the  routine  fixed  the  slope  at  the  current 
value.  However,  at  the  next  time  step  the  moment  exceeded 
the  value  of  MoFtoss,  so  the  slope  was  freed  and  damping 
continued.  The  loop  starts  at  the  far  right  and  ends  on  the 
inner  vertical  line.  The  value  of  Horobs  is  not  reached 
after  the  first  two  damping  cycles  for  the  initial 


FIGURE  4,20  Beam  Shape  vs.  Time 
Second  Half -Cycle 


conditions  given,  so  the  moment  vanes  up  and  down  the  line 
Mhich  indicates  a  constant  slope  o-f  -O.OIOIB  rads. 

■  able  4.1  shows  the  percent  o-f  energy  lost  during  the 
■first  half-cycle  and  a  full  cycle  for  various  clamping 
pressures  and  an  initial  applied  load  of  10,U00  N.  Note 
that  the  results  for  the  first  half-cycle  show  an  increase 
in  damping  as  clamping  pressure  increases.  But,  the  numbers 
reach  a  maximum  and  then  decrease  for  the  damping  done  in 
one  cycleas  clamping  pressure  increases.  This  is  better 
illustrated  in  the  plot  of  Kig.  4.22. 


4.22  Percent  Energy  Loss  vs.  Clamping  iTes 
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CHAPTER  5 

CONCLUSIONS 

the  tinite  element  method  is  capable  o+  predicting  the 
amount  o+  energy  lost  due  to  friction  damping  at  the  joint 
of  a  simple  beam  structure.  It  Mas  able  to  calculate 
slopes,  displacements,  velocities,  and  accelerations,  as 
Mell  as  the  frequencies  associated  with  them.  Enough 
information  was  obtained  that  the  actual  shape  of  the  beam 
could  be  plotted.  The  modified  one— di mensi onal  program 
provided  data  that  appeared  accurate  despite  the  inadequacy 
of  the  original  program  to  handle  the  independent  rotation 
degrees  of  freedom  at  the  joint.  The  problem  was  simplified 
by  the  assumptions  made,  but  the  results  have  opened  the 
door  for  continued  investigation  of  the  friction  damping 
phenomenon. 

A  higher  mode  of  vibration  was  seen  on  the  angular 
velocity  versus  time  plot,  but  the  accuracy  for  the 
frequency  is  debatable  because  only  one  element  was  used. 

The  use  of  only  one  element  made  the  beam  stiffer  than  it 
really  was,  therefore  the  frequency  of  the  second  mode  would 
be  higher  than  it  should  be.  This  could  be  corrected  by 
using  two  elements.  Higher  modes  would  be  introduced,  and 
their  accuracy  would  depend  on  the  number  of  elements  used. 

The  errors  encountered  by  not  using  smaller  time  steps 
and  actually  stepping  over  the  values  of  noROBs  and  the 
angular  velocity  when  turning  the  damping  mechanism  on  and 
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0++  did  not  «++ect  the  results  obtained.  Ihe  hysteresis 
loop  shows  that  the  triction  moment  was  exceeded  rather  than 
equaled,  but  at  the  next  time  step  the  value  at  the  joint  is 
the  value  o+  Mcbrobb.  bo  the  time  step  chosen  was  more  than 
adequate  to  give  accurate  results. 

The  displacement,  Wz,  is  seen  to  decrease  as  the 
clamping  pressure  increases.  This  correlates  to  the  percent 
energy  loss  tor  a  hal-f-cycle  and  a  +ull  cycle.  Ihe  percent 
energy  loss  for  the  first  half-cycle  increases  as  the 
clamping  pressure  increases.  However,  it  reaches  a  maximum 
and  begins  to  decrease  as  clamping  pressure  increases  for  a 
full  cycle. 

The  frequency  increases  as  the  clamping  pressure 
increases.  This  can  be  explained  by  the  beam  becoming 
stiffer  and  the  period  decreasing.  In  other  words,  the  beam 
vibrates  faster.  l-or  a  given  clamping  pressure,  the 
frequency  remains  constant  for  any  initial  applied  load. 

The  amplitude  vanes  with  the  load,  but  the  frequency  does 
not . 

The  accuracy  of  the  results  are  unknown  as  no  data  was 
available  to  compare  with.  However,  the  trends  observed 
indicated  that  the  data  obtained  was  a  good  first  estimate 
for  the  amount  of  friction  damping. 
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RECOMMENDATIONS 

This  thesis  has  opened  the  door  for  many  -follow-on 
studies  to  be  accomplished.  As  many  questions  as  were 
answered,  many  more  new  ones  were  asked.  Some  additional 
work  that  could  be  completed  in  this  area  is: 

1.  A  two-dimensional  -finite  element  program  could  be  used 
to  verify  the  results  obtained  by  the  modified 

one-di mensi onal  program. 

2.  An  investigation  into  the  higher  modes  could  be  made 
using  more  elements  in  the  beam  model. 

3.  The  assumptions  could  be  relaxed  to  allow  more  variables 
to  come  into  play,  such  asj 

a.  Allow  the  clamping  pressure  to  vary  from  a  maximum  at 
the  point  of  application  to  some  lesser  away  from  the  point 
of  application. 

b.  Allow  partial  slip.  This  would  mean  the  joint  would 
always  be  in  some  form  of  slip  once  a  break-away  value  was 
exceeded. 

c.  Allow  the  coefficient  of  friction  to  vary  with  time 


I 


to  simulat*  Msaring  04  tha  contact  sur4acas. 

I 

4.  Finally,  the  problem  should  be  accomplished  ! 

experimentally  to  obtain  some  laboratory  values  to  compare 

against  the  numerical  integration  methods.  C'lvN' 


APPENDIX  A 

Beam  is  made  of  steel 

IX  =  0.3  7827.1093  Kg/m-^^  E  =  2.0E11  N/m^ 

Area  of  cross-section  =  bh  =  (0.012/m) (0.02b4m) 

-  3.22b8E-04  m^ 

^A  *  2.5248689  Kg/m 

1-  <l/12)bh=»  «  (1/12)  (0.0127m)  (0.02b4m)^  -  1 . 7342976E-08  m 
El  «  3468.5952  Nm= 

Area  of  contact  surface  ■  *  H(.012/m>^  ■  5.067074t-04  m 

Mcanosa  «  2n«PK'*/3  »  <1.2b/037t-06  m^)  (P  N/m^)  Nm 


APPENDIX  B 


the  di+ferential  equation  for  the  dynamic  problem  consists 
o+  the  following  terms: 


Mm  Ml-* 
M21  M23  M23  Ma^ 
M3 1  M33  M33  M3^ 
M4 1  Mna  M«3  M.*^ 


Kll  K.13  Ki3  K.i« 
K21  K22  K23  Ksa 
K31  K32  K33  Ks* 
K^»  K-*2  K.*3  K.*a 


The  clamped  boundary  conditions  at  the  -first  node 


(Wi 


Ui  =  u)  allow  the  following  simplification  to  be  made 


M33  M3-* 
M*- 


!  K33  K3.*  1 
K.*.3  k.* 

Therefore,  the  moment  required  to  compare  against  McaFtoss  is 
Ma  ®  M.43W2  -f  M.*«02  -t  K4.3W2  +  K*A02 


The  displacement  and  acceleration  vectors  are  obtained  by 
numerical  analysis,  the  Newmark  Beta  method  in  this  case. 


APPENDIX  C 


SUBROUTINE  STRESS ( NPE , NDF , IBEAM , I ELEM , WO , ELX , N , NEM , I FLAG , NT , 
*VBF, INITIAL,BXO,KICK,FORCE,R,FRICTION, ITEM,IBF,VBDY, IBDY, 
*K0UNT,GF,GF1 ,GF2,T) 


X  . GLOBAL  COORDINATE 

XI  . LOCAL  COORDINATE 

SF . ELEMENT  INTERPOLATION  FUNCTIONS 

GDSF . FIRST  DERIVATIVE  OF  SF  W.R.T.  GLOBAL  COORDINATE 

GDDSF. .. .SECOND  DERIVATIVE  OF  SF  W.R.T.  GLOBAL  COORDINATE 

WO . COLUMN  OF  GENERALIZED  DISPLACEMENTS 

W . INTERPOLATED  GENERALIZED  DISPLACEMENT 

DW . FIRST  DERIVATIVE  OF  W:  DW/DX 

DDW . SECOND  DERIVATIVE  OF  W:  D(DW)/DX 


NOTE:  W,  DW,  AND  DDW  ARE  COMPUTED  AT  NINE  POINTS  OF  EACH 
ELEMENT  ( DW  AND  DDW  ARE  NOT  EXPECTED  TO  BE  ACCURATE 
AT  THE  NODAL  POINTS  OF  THE  ELEMENT) 


IMPLICIT  REAL*8(A-H,0-Z ) 

COMMON/STRS/XG( 1 8 ) ,W( 1 8 ) ,DW( 1 8 ) ,DDW( 18) ,NBDY,NBF, B43 , B44 , 
*G43,G44,VT,VT1 

COMMON/SHP/SF ( 4 ) , GDSF ( 4 ) , GDDSF ( 4 ) , GJ 

DIMENSION  GAUSSO)  ,W0(4)  ,  ELX(  4  )  ,  VBF  (  1  1  )  ,IBF(  1  1  )  ,VBDY(  1  1  ) 
*IBDY( 1 1 ) ,GF ( 4) ,GF1 (4) ,GF2(4) 

DATA  GAUSS/- 1 . ODO , -0 . 7 5D0 ,-0 . 50D0 , 

*-0 . 25D0 , 0  .  ODO  ,  0 . 25D0 , 0 . 50D0 , 

*  0.75D0, 1 .ODO/ 


NET=NPE 

IFdELEM  .EQ.  0)NET=4 
H  =  ELX(NPE)-ELX( 1 ) 

DO  85  NI=1 ,9 
XI  =  GAUSS(NI) 

CALL  SHAPE(XI,H,NPE,NET,IELEM) 
NN=(N-1 )*9+NI 

XG(NN)  =  0.5*H*( 1 .0+XI)+ELX{ 1 ) 

W(NN)=0.0 

DW(NN)=0.0 

DDW (NN) =0.0 

DO  60  1=1 ,NET 

W(NN)  =  W(NN)  +  SF(I)*W0(I) 

DW ( NN ) =DW ( NN ) +GDSF ( I ) *W0 ( I ) 
IFdELEM  .NE.  0)GOTO  60 
DDW ( NN ) =DDW ( NN ) +GDDSF ( I ) *W0 ( I ) 
60  CONTINUE 

IFdBEAM.EQ.O)GOTO  85 
IF(N.NE. 1 .OR.NI.NE.9)GOTO  85 
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